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Abstract 

Neutrinos are the best candidates to test the extreme Universe and ideas beyond the 
Standard Model of particle Physics. Once produced, neutrinos do not suffer any kind of 
attenuation by intervening radiation fields like the Cosmic Microwave Background and 
are not affected by magnetic fields. In this sense neutrinos are useful messengers from the 
far and young Universe. In the present paper we will discuss a particular class of sources 
of Ultra High Energy Cosmic Rays introduced to explain the possible excess of events 
with energy larger than the Graisen-Zatsepin-Kuzmin cut-off. These sources, collectively 
called top-down, share a common feature: UHE particles are produced in the decay or 
annihilation of superheavy, exotic, particles. As we will review in the present paper, the 
largest fraction of Ultra High Energy particles produced in the top-down scenario are 
neutrinos. The study of these radiation offers us a unique opportunity to test the exotic 
mechanisms of the top-down scenario. 



1. Introduction 

The Ultra High Energy (UHE) neutrino detection is one the most important step 
forward in the Cosmic Ray (CR) Physics. The discovery of neutrinos with energy 
larger than 10 17 eV will start the neutrino astronomy, enabling the observation of 
the most distant and powerful sources in the Universe. The existence of neutrinos 
with such high energy is intimately related to the observation of Ultra High Energy 
Cosmic Rays (UHECR). 

Soon after the discovery of the Cosmic Microwave Background (CMB) radiation 
it was shown that the flux of UHECR, with an energy larger than 10 18 eV, should 
be characterized by a sharp steepening at energies ~ 10 20 eV, due to the absorption 
processes on the CMB radiation. This effect is the well known Graisen-Zatsepin- 
Kuzmin (GZK) cut-off [1]. After a few decades of observations the detection of 
the GZK steepening is still one of the major open problems in UHECR physics 
and the experimental data are not conclusive. The 11 Akeno Grand Air Shower 
Array (AGASA) events with energy larger than 10 20 eV [2] contradict the expected 
suppression of the UHECR spectrum. On the other hand the HiRes data seem to be 
consistent with the GZK cut-off picture [2] . If the UHECR primaries are protons 
and if they propagate rectilinearly, as the claimed correlation with BL-Lacs at 
energy 4 — 8 x 10 19 eV implies [3], than their sources must be seen in the direction 
of the highest energies events with energies up to 2 — 3 x 10 20 eV detected by HiRes, 
Fly's Eye and AGASA [2]. At these energies the proton attenuation length is only 
about 20 — 30 Mpc and no counterparts in any frequency band was observed in 
the direction of these UHECR events. This is a strong indication that CR particles 
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with energies larger than 10 20 eV may have a different origin from those with lower 
energies. 

Models of origin of UHECR fall into two categories, top-down and bottom-up. 
In the bottom-up scenario UHECR originates from cosmic accelerators. In these 
accelerators particles of relatively low energy are brought to UHE through multiple 
interactions at the source. The most promising accelerators known in nature are 
based on the diffusive shock acceleration mechanism, in which the particle acceler- 
ation is realized through multiple interactions with a shock front. This mechanism 
works fairly well in Super Nova Remnants (SNR) that are believed to be the re- 
sponsible for the acceleration of Galactic Cosmic Rays of energy E < 10 18 eV (for a 
review see [4]). At the highest energies different bottom- up scenario have been pro- 
posed, among them, the most promising, are those in which acceleration is realized 
through the interaction with a relativistic shock front in Active Galactic Nuclei 
(AGN) and Gamma Ray Bursts (GRB). In the framework of bottom-up models 
the observed UHECR flux should show the predicted GZK steepening and UHE 
neutrinos are produced by the interaction of UHECR with different backgrounds, 
at the source and during their journey to us. These are the so-called cosmogenic 
neutrinos (first proposed in [5]) that we will not discuss here. 

The presence of an excess of events as claimed by AGASA inspired the intro- 
duction of several exotic models for the production of UHECR. These models, 
collectively called top-down, explain the excess of AGASA and give also a clear 
explanation for the lacking of any counterpart of the highest energy events. Many 
different ideas have been proposed among top-down models: strongly interacting 
neutrinos [6] and new light hadrons [7] as unabsorbed signal carriers, Z-bursts [8], 
Lorentz-invariance violation [9], Topological Defects (TD) (sec [10] for a review), 
and Superheavy Dark Matter (SHDM) (see [11] for a review). 

In the present paper we will concentrate our attention on the two last models 
that show common features: UHE particles are produced in the decay of superheavy 
particles, that we shall call collectively X particles, with a typical mass of the order 
of the Grand Unified energy scale Mqut — 10 24 eV. In the case of TD the X particle 
once produced, by the internal dynamics of the defect or through the interaction 
of different defects, immediately decays. While in the case of SHDM the X particle 
itself is long-lived contributing to the Dark Matter of the universe. 

From the point of view of elementary particle physics the X particle decay process 
proceed in a way similar to e + e~ annihilation into hadrons: two or more off-mass- 
shell quarks and gluons are produced and they initiate a QCD cascade. Finally the 
partons are hadronized at the confinement radius. Most of the hadrons in the final 
state are pions and thus the typical prediction of all these models is the dominance 
of neutrinos and photons at the highest energies E > 5 x 10 19 eV. It is important 
to stress here that these models predict neutrino fluxes most likely within reach of 
the first generation neutrino telescopes such as AMANDA, and certainly detectable 
by future kilometer-scale neutrino observatories [12]. 

2. Hadrons spectrum in X decay 

The first step to determine the neutrino flux produced in the decay of X par- 
ticles is the determination of the hadron spectrum. Moreover, this evaluation is 
particularly important because it represents a direct signature of the production 
mechanism that, in principle, can be detected experimentally. As discussed in the 
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introduction, the mass of the decaying particle, Mx, that represents the total CMS 
energy ^/s, is in the range 10 13 - 10 16 GeV. 

The existing QCD Monte Carlo (MC) codes become numerically unstable at 
much smaller energies, e.g., at yfs ~ 10 GeV and the computing time increases 
rapidly going to larger energies. In this section we will briefly review the main results 
obtained, in the computation of the top-down spectrum of UHE particles, using two 
different computational techniques: one based on a new MC scheme [13,14] and the 
other based on the Dokshitzcr-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution 
equations [13,16]. In both cases SUSY is included in the computation. 

Monte Carlo simulations are the most physical approach for high energy calcula- 
tions which allow to incorporate many important physical features as the presence 
of SUSY partons in the cascade and coherent branching. The perturbative part of 
a QCD Monte Carlo simulation is quite standard with the inclusion of SUSY. For 
the non-perturbative hadronization part an original phenomenological approach is 
used in Ref. [13]. The fragmentation of a parton i into an hadron h is expressed 
through perturbative fragmentation function of partons D\ (x, Mx), that represents 
the probability of fragmentation of a parton i into a parton j with momentum frac- 
tion x = 2p/Mx, convoluted with the hadronization functions fj(x,Qo) at scale 
Qo, that is understood as the fragmentation function of the parton i into the hadron 
h at the hadronization scale Qo — 1-4 GeV [13]. To obtain the fragmentation func- 
tions of hadrons one has: 

D?(x,Mx)= ]T f ^D{(- M x )tf(z,Qo) (1) 

Jx z z 

3=1,9 

where the hadronization functions do not depend on the scale Mx- This impor- 
tant property of hadronization functions allows the determination of f^{x, Qo) from 
the available LEP data, D^(x, Mx) at the scale Mx = Mz- 

The fragmentation functions D^(x,Mx) at a high scale Mx can be calculated 
also evolving them from a low scale, e.g. Mx = Mz, where they are known experi- 
mentally or with great accuracy using the MC scheme. This evolution is described 
by the Dokshitzcr-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation [15] which 
can be written as 

dfD H = J- ^M Pij {z) ® D^x/z, t) , (2) 

3 

where t = ln(s/s ) is the scale, <g> denotes the convolution f®g = J z dx/xf(x)g(x/z), 
and Pij is the splitting function which describes the emission of parton j by parton 
i. Apart from the experimentally rather well determined quark fragmentation func- 
tion Dg(x, Mz), also the gluon fragmentation function Dg(x, Mz) is needed for 
the evolution of Eq. (2). The gluon FF can be taken either from MC simulations 
or from fits to experimental data, in particular to the longitudinal polarized e + e~ 
annihilation cross-section and three-jet events. The first application of the DGLAP 
method for the calculation of hadron spectra from decaying superheavy particles 
has been made in Refs. [16]. The most detailed calculations have been performed 
by Barbot and Drees [16], where more than 30 different particles were allowed to 
be cascading and the mass spectrum of the SUSY particles was taken into account. 
The results obtained with the two different techniques discussed above agree fairly 
well [13]. The accuracy in the hadron spectrum calculations has reached such a 
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level that one can consider the spectral shape as a signature of the model. The 
predicted hadron spectrum is approximately oc dE/E 1 - 9 in the region of x relevant 
for UHECR observations. 



3. Spectra of Neutrinos, Photons and Nucleons 

The spectra of neutrinos and photons produced by the decay of superheavy parti- 
cles are of practical interest in high energy astrophysics and can be computed from 
the decay of charged pions [13]. The FFs for charged pions and protons+antiprotons 
can be determined, following [13], from the FFs of hadrons Dh simply introducing 
the ratios En(x) and e 7T (x) as: Djy(x) — SN{x)Dh{x) and D 7T (x) = e v (x)Dh(x). 
The spectra of pions and nucleons at large Mx have approximately the same shape 
as the hadron spectra, and one can use in this case — 0.73 ± 0.03 and Em = 
0.12 ±0.02 [13]. 

An interesting feature of the up-dated calculations performed in [13] and by 
Barbot and Drees in [16] is the ratio of photons to nucleons, j/N. At x ~ 1 x 10~ 3 
this ratio is characterized by a value of 2 - 3 only [13]. This result has an important 
impact for SHDM and topological defect models because the fraction of nucleons 
in the primary radiation increases. However, in both models photons dominate (i.e. 
their fraction becomes > 50%) at E > (7 - 8) x 10 19 eV. 

Let us now concentrate our attention on UHECR from superheavy dark matter 
(SHDM) [18] and topological defects (TD) [19]. The comparison of the UHECR 
spectrum obtained with the AGASA data, will provide us with the correct neutrino 
flux normalization. 

Production of SHDM particles naturally occurs in a time-varying gravitational 
field of the expanding universe at the post-inflationary stage. The relic density of 
these particles is mainly determined (at fixed reheating temperature and inflaton 
mass) by their mass Mx- The range of practical interest is (3 — 10) x 10 13 GeV, 
at larger masses the SHDM is a subdominant component of the DM. SHDM is 

—halo —halo 

accumulated in the Galactic halo with the overdensity S = ^&tr = n DM — , where 

Pdm ~ °- 3 GcV/cm 3 , p cr = 1.88 x 10~ 29 /i 2 g/cm 3 and CDM /i 2 = 0.135 [20]. With 
these numbers, 8 w 2.1 x 10 5 . Because of this large local overdensity, UHECRs from 
SHDM have no GZK cutoff. 

dumpiness of SHDM in the halo can provide the observed small-angle clustering. 
The ratio rx = ^x{to/Tx) of relic abundance Ox and lifetime tx of the X particle 
is fixed by the observed UHECR flux as rx ~ 10~ n . In the most interesting case 
of gravitational production of X particles, their present abundance is determined 
by their mass Mx and the reheating temperature Tr. Choosing a specific particle 
physics model one can fix also the life-time of the X particle. There exist many 
models in which SH particles can be quasi-stable with lifetime t x S> 10 10 yr. 
The measurement of the UHECR flux, and thereby of rx, selects from the three- 
dimensional parameter space (M X ,T R , tx) a two-dimensional subspace compatible 
with the SHDM hypothesis. 

In Figure 1 (left panel) we have performed a fit to the AGASA data using the 
photon flux from the SHDM model and the proton flux from uniformly distributed 
astrophysical sources. For the latter we have used the non-evolutionary model of 
[21]. The photon flux is normalized to provide the best fit to the AGASA data at 
B>4x 10 19 eV. 
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Fig. 1. [Right Panel] Comparison of SHDM prediction with the AGASA data. The calculated 
spectrum of SHDM photons is shown by the label 7 and by the label p e xtr the spectrum of 
extragalactic protons (uniformly distributed astrophysical sources). The sum of these two spectra 
is shown by the thick black curve. The red thick line is the SHDM neutrino flux. [Left Panel] 
Diffuse spectra from necklaces. The red thick curve shows neutrino flux, the black thick curve is 
the sum of protons and photons fluxes produced by necklaces (labeled thin black lines). 



One can see from the fit of Figure 1 (left panel), that the SHDM model can explain 
only the excess of AGASA events at E > 1 X 10 20 eV: depending on the SHDM 
spectrum normalization and the details of the calculations for the extragalactic 
protons, the flux from SHDM decays becomes dominant only above (6—8) x 10 19 eV. 

Topological Defects (for a review see [22] and reference therein) can naturally 
produce UHE particles. The following TD have been discussed as potential sources 
of UHE particles: superconducting strings, ordinary strings, monopolonium (bound 
monopole-antimonopole pair), monopolonia (monopole-antimonopole pairs connected 
by a string), networks of monopoles connected by strings, vortons and necklaces (see 
Ref. [22] for a review and references). Monopolonia and vortons are clustering in 
the Galactic halo and their observational signatures for UHECR are identical to 
SHDM. However the friction of monopolonia in cosmic plasma results in monopolo- 
nium lifetime much shorter than the age of the universe. Of all other TD which are 
not clustering in the Galactic halo, the most favorable for UHECR are necklaces. 

Necklaces are hybrid TD produced in the symmetry breaking pattern G — > H x 
U(l) — > H x Z<i. At the first symmetry breaking monopoles are produced, at the 
second one each (anti-) monopole get attached to two strings. This system resembles 
ordinary cosmic strings with monopoles playing the role of beads. Necklaces exist 
as the long strings and loops. The symmetry breaking scales of the two phase 
transitions, r\ m and 77,,, are the main parameters of the necklaces. They determine 
the monopole mass, m ~ Air'q m /e 1 and the mass of the string per unit length \i ~ 
27TT7g . The evolution of necklaces is governed by the ratio r ~ m/fid, where d is the 
average separation of a monopole and antimonopole along the string. As it is argued 
in Ref. [23], necklaces evolve towards configuration with r > 1. Monopoles and 
antimonopoles trapped in the necklaces inevitably annihilate in the end, producing 

heavy Higgs and gauge bosons (X particles) and then hadrons. The rate of X 

2 

particles production in the universe can be estimated as [23] fix ~ Wm > wnere t 
is the cosmological time. 

The photons and electrons from pion decays initiate e-m cascades and the cascade 
energy density can be calculated as w cas = \f- K r 2 [i J Q ° ^ n+ z )i = I/tt 7 " 2 ^, where 
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z is the redshift and f w ~ 1 is the fraction of the total energy release transferred to 
the cascade. The parameters of the necklace model for UHECR are restricted by 
the EGRET observations [24] of the diffuse gamma-ray flux. This flux is produced 
by UHE electrons and photons from necklaces due to e-m cascades initiated in colli- 
sions with CMB photons. In the range of the EGRET observations, 10 2 — 10 5 MeV, 
the predicted spectrum is oc E~ a with a — 2 [25]. The EGRET observations de- 
termined the spectral index as a = 2.10 ± 0.03 and the energy density of radiation 
as a; bs ~ 4 x 10 -6 eV/cm 3 . The cascade limit consists in the bound tu cas < lu \> s . 
According to the recent calculations, the Galactic contribution of gamma rays to 
the EGRET observations is larger than estimated earlier, and the extragalactic 
gamma-ray spectrum is not described by a power-law with a = 2.1. In this case, 
the limit on the cascade radiation with a = 2 is more restrictive and is given by 
w C as — ^ x 10 _6 eV/cm 3 ; we shall use this limit in further estimates. Using w cas with 
f n = 1 and to = 13.7 Gyr [20] we obtain from the limit on the cascade radiation 
r 2 fi < 8.9 x 10 27 GeV 2 . 

The important and unique feature of necklaces is their small separation D, which 
ensures an high density. The distance D is given by D ~ r -1 / 2 ^ [23]; since r 2 fi 
is limited by e-m cascade radiation we can obtain a lower limit on the separation 

between necklaces as D ~ ( j^f h > 10(>/10 6 GeV 2 ) 1/4 kpc, this small 
distance is a unique property of necklaces allowing the unabsorbed arrival of parti- 
cles with the highest energies. The fluxes of UHECR from necklaces are shown in 
Figure 1 (right panel). We used in the calculations r 2 /i = 4.7 x 10 27 GeV 2 which 
corresponds to w cas = 1.1 x 10 -6 eV/cm 3 , i.e. twice less than allowed by the bound 
on uj cas . The mass of the X particles produced by monopolc-antimonopolc annihi- 
lations is taken as Mx = 1 x 10 14 GeV. From Figure 1 (right panel) one can see 
that the necklaces model for UHECR can explain only the highest energy part of 
the spectrum, with the AGASA excess somewhat above the prediction. Thus UHE 
particles from necklaces can serve only as an additional component in the observed 
UHECR flux. This result has a particular impact on the possible UHE neutrino 
detection. In fact, the necklaces model is only under constrained by the available 
UHECR data, in this context only a clear UHE neutrino observation with a typical 
spectrum as in figure 1 (right panel) can confirm (or falsify) the model. 
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